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Abstract 

A quantum control landscape is denned as the observable as a function(al) of the system con- 
trol variables. Such landscapes were introduced to provide a basis to understand the increasing 
number of successful experiments controlling quantum dynamics phenomena. This paper ex- 
tends the concept to encompass the broader context of the environment having an influence. 
For the case that the open system dynamics are fully controllable, it is shown that the control 
landscape for open systems can be lifted to the analysis of an equivalent auxiliary landscape of 
a closed composite system that contains the environmental interactions. This inherent connec- 
tion can be analyzed to provide relevant information about the topology of the original open 
system landscape. Application to the optimization of an observable expectation value reveals 
the same landscape simplicity observed in former studies on closed systems. In particular, no 
false sub-optimal traps exist in the system control landscape when seeking to optimize an ob- 
servable, even in the presence of complex environments. Moreover, a quantitative study of the 
control landscape of a system interacting with a thermal environment shows that the enhanced 
controllability attainable with open dynamics significantly broadens the range of the achievable 
observable values over the control landscape. 



*Departmcnt of Chemistry, Princeton University, Princeton, New Jersey 08544, USA 



1 



I. INTRODUCTION 



Quantum optimal control concepts are being used to address many applications in- 
cluding bond-selective chemical reactions, quantum state engineering, and quantum com- 
putation , |6j. Toward the natural goal of finding the best attainable outcome, 
optimal control theory can be applied to redirect the quantum dynamics subject to any 



3]. In the laboratory, optimal control 
8[ that guide a closed- 



practical constraints on the laboratory controls 
experiments generally have been based on learning algorithms 
loop search for an effective control despite the lack of full information about the system's 
Hamiltonian, its interaction with the environment and the influence of noise. The growing 



number of experimental successes 



ll| collectively show that it is (unexpectedly) 



easy to find high-quality control results even with severely constrained control fields (i.e., 
shaped pulses from the same Tiisappire laser operating at ~800nm with a bandwidth of 
~20nm are used as the driver for most of the experiments). 

The quantum control landscape is the expectation value of the physical observable as 
a functional J = J[e] of the time dependent control field e(-). Optimal controls maximize 
(or minimize, depending on the control goal) the objective functional, i.e., such that 
J[e] = max e J[e] (resp., J[e] = min e J[e]). In the laboratory, for controlling a complex 
system typically detailed knowledge about the dependence of the objective functional on 
the control is not available, and the practical search for an optimal control solution is 
performed adaptively via various algorithms (e.g., genetic algorithms, gradient methods). 
If the objective functional has local traps (local minima or maxima), then their presence 
can slow down the efficiency of the search or even permanently trap the algorithm in 
a local maximum (resp., minimum). This circumstance motivates the analysis of the 
topological properties of quantum control landscapes. In particular, a priori information 
about the absence of local traps for common quantum objective functionals for a wide 
class of quantum systems would imply that the laboratory searches would not be hindered 
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and a global optimal solution would be finally reached (assuming that the control field is 
sufficiently flexible) . Moreover, the lack of landscape traps would provide an explanation 



'or the evident re 
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ative practical ease of obtaining optimal solutions. Recent studies 



13, 
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191 ] show that such favorable landscape properties exist for 



optimizing the expectation value of a system observable for any controllale finite level 
closed quantum system. 

Realistic quantum systems are always exposed to some kind of environment (e.g., 
chemical reactions in heat baths, or atoms in optical cavities), and their landscape analy- 
sis for such systems is important. In most cases, the environmentally induced irreversible 



quantum dynamics 
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21] would downgrade the quality of the control outcome, espe- 



cially when little can be done to tailor the environment. This situation can be serious in 
many applications, especially quantum information sciences [2l|]. However, there are in- 
dications that it may be possible to fight against or even cooperate with environmentally 



induced dynamics 



22] . Moreover, combining incoherent control by the environment with 



coherent control by an electromagnetic field offers a general technique for manipulating 
quantum systems by optimally affecting both the Hamiltonian and dissipative aspects 
of the dynamics; the latter circumstance corresponds to the environment consisting of 
incoherent radiation or a gas of electrons, atoms or molecules whose state is suitably op- 
timized H (e.g., to allow for the creation of mixed target states). There also have been 
explicit applications of controlled quantum correlations induced to non-classically ma- 
nipulate quantum states (e.g., through quantum error correction using redundant qubits 
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25(1, atomic control in single- mode optical cavities [26|], and real-time feedback con- 
trol 27], |28j). These scenarios show the prospects of attaining full control of the system 
evolution by "freezing" orpossibly counter-intuitively "reversing" what was thought to 



be irreversible dynamics 



29J. In the laboratory, even if an explicit treatment of the 



environment is not considered, under optimization of the observable the control will nat- 
urally address the environment to maximally draw on its beneficial features or minimize 
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its deleterious effects to best achieve control of the system. 

This perspective forms the background for the present paper aiming to investigate con- 
trol landscapes of open quantum systems whose dynamics can be effectively manipulated. 
The landscape analyses will be explored within the framework of the Kraus operator-sum 



representation of open system dynamics 



301 ]. revealing that no local suboptima exist in 



the control landscape of open systems, although there is an increase of the number of sad- 
dle submanifolds from that in the case of analogous closed quantum systems. The open 
system landscape analysis indicates that the search for optimal controls in the laboratory 
may not be significantly hindered by the presence of an environment, thereby providing 
a basis to understand the relative ease of obtaining optimal control of open systems and 
allowing for the development of more effective searching algorithms. 

The paper is organized as follows. Section II summarizes the properties of Kraus maps, 
while section III presents the concept of landscape lifting and discusses its general aspects. 
Section IV applies landscape lifting to Kraus maps for manipulation of an observable, 
thereby revealing a rich structure for the critical submanifolds in the search space and 
the quantitative influence of the environment on the landscape. Section V summarizes 
the results. 



II. REPRESENTATIONS OF OPEN QUANTUM SYSTEM DYNAMICS 

We summarize the Kraus operator-sum formalism of open quantum system dynam- 
ics via a simple "system plus environment" model, where the quantum system with 
Hilbert space dimension N is coupled to an external quantum environment (taken as A- 
dimensional) , and the composite of the system and the environment obeys the Schrodinger 
equation, 

ih ^° tal = [-Htotal, Ptotal], (1) 
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where ptotai is the total density matrix. The total Hamiltonian H tota x includes the in- 
ternal Hamiltonians of the system and the environment as well as their interaction, all 
contributing to the total system evolution operator U tota \(t) on the total Hilbert space 
TC = Hs <S>TCe, where TCs and He are the Hilbert spaces of the system and environment, 
respectively. The composite system is assumed to initially be prepared in a product 
state ptotai(O) = p <8> g, where p and g are the initial states of the system and the en- 
vironment, respectively. The system dynamics described by ps{t) can be obtained by 
tracing p to tal(*) = f / totai(^)Ptotai(0)L 7 t t otal (t) over the environment, leading to the Kraus 
operator-sum representation: 

A 

p s (t) = K ^(t)pK^(t), (2) 

a,p=l 

where K a p{t) =Ti E {U tota i(t)(I N ® g 1/2 \f3)(a\)} is the (a, (3)-th NxN block of the XN x 
XN matrix K = U to t a ,i(t)(lN®g 1 ^ 2 ) under some arbitrary orthonormal basis {\a), • • • } 
of H.E- The matrices {K a p} form the Kraus representation of the dynamical map ([2]). 
For compactness, we use K to denote the Kraus representation by the set of Kraus 
operators {K a p} associated with the initial environmental state g, which always satisfies 
the identity 

K ] K = I N ®g (3) 

to guarantee that the trace-preservation of the system density matrix jsij ]. The set of 
Kraus representations induced by the environment state g is denoted by K[g, N], which 
can be taken as the image of the unitary group IA(\N) under the mapping 

K = & e {U) = U(I N ®g l/2 ), UeU(XN). (4) 

Suppose that the rank of g is p, then 1K[^, N] is homeomorphic to the following homoge- 
neous space of U(XN) 



where each K e K[^, N] corresponds to a coset in the unitary group. 

The definition ([3]) of Kraus representation is still ambiguous because the reduction (T2J) 
to the system density matrix is dependent on the choice of basis for He- In general, two 
Kraus operators K and K correspond to the same Kraus map if 

a,/3 a 1 , (3' 

It is easy to show that the following set of Kraus representations 

{K = (I N ® V)K, V eU{\)} (6) 

are equivalent to K, and they exhaust all equivalent representations when q is a pure 
state. Here the unitary matrix V represents all possible basis transformations on the 
environmental subsystem. 

Since the sequential order of the Kraus operators in the sum does not affect the opera- 
tion of Kraus maps on the density matrix, and the effective number of independent Kraus 
operators is at most N 2 because the Kraus maps are linear in the space of iV-dimensional 
Hermitian operators, any Kraus map can be associated with an iV 2 -dimensional "ficti- 
tious" environment starting from a pure state. This characterizes the whole set of open 
system Kraus maps without any explicit environmental details appeared in the formu- 
lation. In more general cases, ~K[g, N] has to be associated with a specific environment 
initial state g and the set only involves a proper subset of all possible dynamical maps of 
the underlying open quantum system. 

III. LANDSCAPE LIFTING 

The quantum control landscape is defined by a cost functional </[£(£)] built on the space 
of time-dependent control fields that enter into the Hamiltonian H tota \ to manipulate the 
quantum dynamics. A main objective of the analysis is to determine the topology of the 
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landscape extremal controls (i.e., critical points of the landscape function) and understand 
the influence of the topology upon the search for optimal controls. 

It can be difficult to carry out a landscape analysis in the space of time-dependent 
control fields that form a subset of the infinite dimensional space £ 2 (M) of Lebesque 
integrable functions. However, under some special circumstances with symmetries, the 
arguments of the functional J(-) may be expressed as a function of a smaller group of 
variables to define the landscape. For example, consider the state-to-state quantum tran- 
sition landscape J[e(-)] = | (f\U[e(-), tf\ \i) | 2 , where U[e(-),tf] is the system propagator 
at time tf under control e(-) at going from the initial state \i) to the target state |/). 
The landscape function can be rewritten as J{U) = \{f\U\i)\ 2 over the iV 2 -dimensional 
unitary transformation group manifold U(N), or more compactly, J{ip) = \(f\ip)\ 2 over 
the unit sphere in the iV-dimensional Hilbert space of quantum states \ip). A general 
connection between such landscapes is shown by the following theorem: 

Theorem 1 Let L(x) be a landscape defined on some differential manifold X. If x is a 
function x = h(y) of y in some other differential manifold y with dim y > dim X, and 
h is locally surjective on a neighborhood of some point y e y , i.e., the Jacobian has full 
rank: 

d h 

= dimX 



, dh 
rank - — 
dy 



(7) 

x=h(y ) 



y=yo 

in some local coordinate system. Theny is critical for Loh in y if and only if x Q = h(y ) 
is critical for L in X , and they have an identical number of positive and negative Hessian 
eigenvalues at y and x , respectively. Moreover, if the inverse image h~ l (x Q ) of every 
critical point x is connected, then the connected components of their critical manifolds 
are one-to-one between the two landscapes. 



Proof: From the chain rule 

d(Loh) _ dL dh 
dy dx dy' 



(8) 
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the rank condition (J7j) guarantees that the gradient of L o h vanishes at y if and only if 
the gradient of L vanishes at Xq, hence yo is critical if and only if xq is critical. Moreover, 
the Hessian quadratic forms at the critical points are related by: 

d 2 {Loh) _ /d/A f d 2 L dh 
dy 2 \dy J dx 2 dy 

The right hand side is a congruent transformation on the Hessian quadratic form of L 
at x under the same full-rank condition ([7]), which preserves the signs of the nonzero 
Hessian eigenvalues. Provided that the Hessian form is sufficient to determine the types 
of critical points for both functions, then any maximum (minimum, saddle) point in y 
will be mapped to a maximum (minimum, saddle) point in X. 

Every connected critical manifold of L o h in y must produce a connected one of L in 
X because the mapping h : X — > y is continuous, so the number of critical components 
in X must be no larger than those in y. In a general setting, the number of connected 
components and their simple (multiple) connectedness may not be preserved. Moreover, 
let 2/1 and y 2 belong to two connected critical submanifolds in y, respectively, that are 
mapped into the same point xq in a connected critical component in X. They must lie 
within the connected component h~ 1 (xo), which is connected by assumption. Therefore, 
the two critical submanifolds in y can be connected, and hence every connected compo- 
nent in X must be mapped from a unique connected critical submanifold in y, i.e., the 
connected critical submanifolds in X and y are one-to-one. End of proof. 

The rank condition is essential to the viability of lifting the landscape analysis results 
in a kinematic picture (i.e., on the unitary group irrespective to the system dynamics) 
to the control landscape on the space of control fields, which has been adopted in prior 
studies 12, [k], Q. Here, the full-rank property corresponds to the local controllability 
along the trajectory driven by some reference control field. If the reference control steers 
the system propagator to some (kinematic) landscape critical point in U(N) and local 
controllability holds, it must also be critical and possess the same optimality status. 

I L 

Such controls are conventionally called regular extremal controls [32]). According to ([51 , 
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a rank-deficient landscape lifting can lead to a critical control field as well, which is called 
a singular extremal control that may have no correspondence to any kinematic critical 
point in U(N). So, the critical topology in the space of control fields depends on both 
regular and singular extremal controls. The optimality of singular extremals can generally 



be excluded 



32 



331 ]. and hence does not affect the search of optimal controls. Therefore, 



it is sufficient to consider only the normal extremal controls as the inverse image of 
the kinematic critical points in U(N) via landscape lifting. This greatly simplifies the 
landscape analysis by merely working in the kinematic picture. 



IV. QUANTUM ENSEMBLE LANDSCAPES OF OBSERVABLE PREPARA- 
TION 

Theorem [TJ provides a useful means to obtain the critical submanifold of an unknown 
landscape by identifying it with the image of some known landscape via a faithful lifting 
that satisfies the conditions in Theorem 1. Consider the goal of maximizing the expecta- 
tion value of the observable 9 of an open quantum system with initial state p, which can 
be formulated on the set of Kraus operators: 

J(K) = Tr ( K^pKlA = Tr{K(p ® I x )tf{9 ® J A )}. (9) 

\a,/3=l / 

We can lift this landscape via the mapping to the following landscape over IA(\N): 

J(U) = Tr {^ e (U)(p <g> h)^(U)e} = Tr (UPU^Q) , (10) 

where P = p® g and = 9 (g> I\. Because K[g, N] is homeomorphic to the homogeneous 
space of U(XN) as shown in (jSj), the properties in the theorem above are naturally guar- 
anteed from the facts that (a) the mapping of a canonical projection is locally surjective 
3^ | and (b) the inverse image of every Kraus operator K is homeomorphic to its coset in 
U(XN) and hence is connected. Thus, complete information can be extracted about the 
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critical landscape topology for open systems from the auxiliary landscape ffTUl) . which en- 
compasses the composite "system plus environment" model that targets the maximization 
of the expectation value of the same system observable (minimization may be treated just 
as well). We assume that the system is controllable over the set of Kraus maps K[p, N] 
generated from an environmental initial state g, while the controllability over U(XN) is 
not necessarily as demanding as it might seem in practical physical circumstances since 
only a portion of the composite system (e.g., just the system and its immediate local 
environment) is likely required for manipulation. Furthermore, most realistic laboratory 
control goals will admit as fully acceptable just reaching some reasonable neighborhood 
of the absolution maximum objective yield. 

Note that the controllability criterion does not give any information about the control 
landscape for open quantum systems, which might be extremely complex in presence of 
environmental interactions. However, an inherent property revealed about the landscape 
for open systems can be immediately observed from the above landscape lifting analysis. 
The open system landscape is always devoid of false traps whenever the open system 
dynamics is fully controllable, because the equivalent landscape for the closed composite 
quantum system was proven to always possess this basic novelty [if], [iij owing to the 
linearity and unitarity of quantum system dynamics. This result is of fundamental im- 
portance for understanding the evident ease of identifying effective control fields in the 



presence of uncertainties 



10 



11 



12|, Q, Q through the guarantee of not being 



caught by any false traps using suitable search algorithms in the laboratory with flex- 
ible control fields. In addition, the ability to decrease the entropy in the system state 
shows that open system controlled dynamics can outperform control of the corresponding 
isolated system without introducing any false traps, as will be seen later. 

The details of the landscape topology is affected by the initial environmental state, 
which can be resolved via the equivalent landscape (Tl0|) as well. Other work [2] developed 
the technique of contingency tables to characterize critical submanifolds. To be specific, 
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suppose that P has r distinct eigenvalues with degeneracy degrees mi, • • • , m r and 9 has s 
distinct eigenvalues with degeneracy degrees n\ , ■ ■ ■ ,n s . Then every critical submanifold 
of (fTUj) corresponds to a contingency table shown in Table I. 





rtl ■ 


■ n s 


mi 


Cll ' 




m r 


Crl • 


c rs 



Table I: The contingency table of the composite system. The entries satisfy the constraints 
J2i Cij = rrij and ]T\ <Hj = n*. 

Hence, the determination of critical submanifolds can be reduced to a combinatorial 
problem of enumerating all possible contingency tables. The dimension and the numbers 
of positive and negative eigenvalues of the Hessian matrix can also be computed by these 
contingency tables. 

Considering the composite system plus an environment, the marginal conditions can be 
determined by the degeneracy structures of p, g and 9. Suppose that p has r distinct non- 
negative eigenvalues 1 > pi > ■ ■ ■ > p r -i > p r = with multiplicities d\, • • • ,d r (d r can be 
zero corresponding to no zero eigenvalues in p); g has m eigenvalues 1 > qi > ■ • ■ > q m — 
with multiplicities /i, • • • , f m (fm can be zero corresponding to no zero eigenvalues in 
g). The distinct eigenvalues of 9 are r% > • • • > r s with multiplicities ex, ■ ■ ■ ,e s . For 
simplicity, we assume no accidental degeneracies in p <g> g, i.e., piqj ^ pkqe for arbitrary 
1 < i 7^ k < r or 1 < j 7^ i < m. In this case, the degeneracy structure of the initial 
state P = p <S> g is 

Dxi = dxfx, ■ • • , Dij = difj • ■ • , _D r _i im _i = f^-i/m-ii 
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corresponding to nonzero eigenvalues p^-, where i = 1, • • • , r and j = 1, • • • , m, while 
the zero eigenvalue of P has multiplicity D rm = d r \ + f m N — d r f m . The observable 
= 6 ® J A possesses the same group of eigenvalues as 6* whose degeneracy indices are 
multiplied by A, i.e., 

Ei = Aei, • • • , E s _i = Ae s „ 1? E r = \e s . 

The degeneracy structures define the contingency Table II, where the row and column 
sums to be satisfied are (Du, • • • , -D r _i im _i, D rm ) and (E x , • • • , E s ). 





Ei 


£^2 




D u 


cm 


C112 


Clls 


Dr—l,m—l 


Cj — l,m— 1,1 


Cr-l,m-l,2 • 


Cj — l,m— l,s 


D rm 


Crml 


Crm2 


Crms 



Table II: The contingency table of the composite system. 



A. The Full Control Landscape of Kraus Maps 

Consider the control landscape of all possible Kraus maps. As discussed previously, 
this situation can be equivalently associated with a special A 2 -dimensional environment 
which is initially in a pure state g that has only one nonzero eigenvalue. The associated 
contingency Table III shows that P has the same group of eigenvalues as p, except for 
the increase in the number of zero eigenvalues by N 3 — N. 
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N 2 ei ■ 


• N\ s 


d 1 




en • 


■ C U 


d r -i 




Cr-1,1 • 


Cr— l,s 


d r + {N 3 - 


N) 


Crl • 


C rs 



Table III: The contingency table for the landscape with pure environmental initial states. 

Because of the marginal conditions, the first r — 1 rows uniquely determine the con- 
tingency table provided that the sum of the first r — 1 elements in each column does not 
exceed the corresponding column sum. This criterion is automatically satisfied in Table 
III because 

Cij + h Cr-ij < di + h d r -i < N < N 2 6j, j = 1, • • • , s. 

Thus, the column-sum constraints are actually ineffective, which implies that the number 
M of contingency tables depends only on the number of ordered partitions of d±, ■ • • , d r -i, 
i.e., 



r-1 



A/~= TT (11) 

The dependence of Af on the degeneracy indices of the nonzero eigenvalues of p manifests 



itself in the contribution of the entropy of the system to the landscape topology. Appar- 
ently, the diverse nature of the initial state distributions tends to give rise to more critical 
submanifolds for large N. The observable 6 affects the number of critical submanifolds 
only by its number of distinct eigenvalues. A physical interpretation of this number is 
the resolution required to measure the target observable, whose low value facilitates the 
control and measurement, and hence reduces the complexity of the control landscape 
despite the environment being present. 

Concerning the local topological landscape properties, a prime interest is the number 
of nonzero eigenvalues at the global optima, which represents the number of effective 
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search directions close in that neighborhood. The associated contingency Table IV has 
nonzero entries only in its first column and last row. Applying the formulas developed in 
[3], we can calculate the number of nonzero Hessian eigenvalues (all negative) as 

M- = 2N 2 (N -dr)(N - ei). 

This characteristic is proportional to the number of zero eigenvalues in p and the number 
of eigenvalues smaller than e\ in the observable 9. Again, the low entropy of p helps 
to reduce the number of search directions towards the global optima, and improves the 
robustness of this critical manifold. Similarly, the large size e\ of the eigenspace of the 
maximal eigenvalue of 9 also facilities the search near the optima. 





N 2 ei 


■ N 2 e s 


dx 






■ 


d r -\ 




d r -i 


• 


d r + (N 3 - 


-N) 


N 2 e X -N + d r ■■ 


• N 2 e s 



Table IV: The contingency table for the global extremum of the open system control landscape. 



B. Control Landscapes with a Thermal Environment 

This section considers the control landscapes under the common situation that the 
environment is initially in a thermal equilibrium state 

where T is the temperature and k is the Boltzmann constant. For simplicity, we suppose 
that the environmental Hamiltonian has a non-degenerate spectrum, i.e., fx — ■ ■ • — 
f m = 1. Although under most circumstances the associated Kraus representations do 
not cover all possible maps, the landscape still exhibits no false traps according to the 
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analysis. Appendix A estimates the approximate numbers of critical submanifolds un- 
der different circumstances. In the limit of zero temperature, the number reaches its 
maximum value as shown in equation (jlip . and remains constant when A is sufficiently 
large. As the temperature increases, more critical submanifolds grow out from the rela- 
tively flat low-temperature landscape. The set of critical submanifolds are identical for 
all finite temperature values and their number increases exponentially with the size of 
the environment, while the geometry (e.g., the curvature near the critical submanifolds) 
varies with the temperature. In the limit of infinite temperature, the landscape critical 
submanifolds merge into a smaller group whose number increases only polynomially with 
the size of the environment. 

Open system control can also extend the scope of manipulation over the expectation 
value of the target observable, which is reflected by the dynamical range of the control 
landscape defined as A J = J max — J m i n . For closed systems that evolve from pure states, 
the attainable landscape values range from the maximal eigenvalue r% of the target ob- 
servable 9 to its minimal eigenvalue r s , i.e., AJ — r± — r s ; the range is usually further 
restricted in the case of mixed states and in the presence of external disturbances 35]. 
The dynamical range of open quantum systems may be enhanced through control over 
the environment and its interactions with the system. Fig.l examines the dynamical 
range of a two-level quantum system, where the observable is chosen as the Pauli oper- 
ator a z whose eigenvalues are r\ = 1/2 and r 2 = —1/2. The landscape always has the 
ideal dynamical range A J max = 1 at zero temperature, and shrinks when the temperature 
increases. Comparing the trends for coupling to different environments, we see that the 
dynamical range is wider when the open dynamics is fully controllable and when interact- 
ing with larger environments. In the limit of infinite temperature, the dynamical range 
approaches the lower bound for that of the isolated system. In that case the enhanced 
controllability of the open dynamics benefits the control landscape very little. The dis- 
tinction between the two groups of curves that differ by the initial system state p shows 
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that lower-entropy helps to extend the dynamical range, i.e., it facilitates higher-quality 
controls. The degree of attainable control over the environment will be an application 
specific matter. Importantly, it is reasonable to expect that only the neighboring envi- 
ronment to the system will likely need to be controlled, and this prospect may open up 
practical conditions for effective use of the environment to aid in the control of system 
dynamics. 
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Fig 1: Dynamical range of a two-level system in different thermal environments including 
one and two harmonic oscillators. The eigenvalues of 6 are r\ = 1/2 and T2 = —1/2. The 
populations of p with higher entropy are p\ = 0.55 and p2 = 0.45, while p\ = 0.70 and p2 = 0.30 
for p with lower entropy. 

We illustrate the behavior of the search for optimal dynamical control of the above 
two- level spin-| system interacting with some relatively small environments represented 
by a spin-^- system (here we choose A = 6). The total system Hamiltonian is H tota _\ = 
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Lower von Neumann entropy 
Higher von Neumann entropy 



H + e(t)a x <g> I x with 

H = UJ a z ® J A + J 2 ® ^ e J 2 +7 (<T a <g> J Q ), 

ce=x,y,z 

where the cr's are the standard Pauli matrices and the spin operators J x , J y , J z form a basis 
of the 50(3) Lie algebra and possess a A-dimensional irreducible unitary representation 
of the 50(3) Lie algebra basis. The environment is initially in a thermal state (|12[) . where 

Although the control field only acts directly on the system itself, it can be demon- 
strated that the rank of the Lie algebra spanned by H and a x <E> Ix equals that of the 
Lie algebra of U(2X), which implies that the composite system is controllable via the 



361 ] . Fig. 2 shows the 



system-environment coupling given a sufficiently large final time 
maximization of the cost function using the conjugate gradient algorithm, and the 
resulting system entropy at each iteration step. As already shown in Fig.l, the dynamical 
search finds fields that transfer population to the target level beyond the capability of the 
closed system control landscape at different temperatures. Moreover, the optimal control 
fields also purify the system as shown in the decrease of the von Neumann entropy in 
Fig. 2. Hence, a clever use of the environment may improve the accessible dynamical range 
on the system landscape. Fig. 3 shows the resulting optimal control field with and with- 
out the presence of the environment (i.e., for the case of T = 300i^). With a relatively 
small environment, the resonant mode for the state transition arising from the system 
free Hamiltonian is still dominant in the optimal control field. The entropy in the sys- 
tem state drops slightly more quickly when the temperature of the environment is higher 
before optimization, i.e., the controlled system becomes cooler when initially bathed in a 
hotter environment. Although such behavior are physically allowed, the phenomenon is 
still counter-intuitive. This may be an algorithm and objective-dependent issue with the 
present goal maximization of the expectation value of 9 without any specific consideration 
of minimizing the system entropy. 
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Fig 2: Dynamical search for optimal controls in a two-level open quantum system, where the 
dimension of the associated environment is A = 6. The observable is 9 = a z . 

V. CONCLUSION 

This work presents an investigation of control landscape topology for open quantum 
systems. The analysis was performed with the aid of landscape lifting and reveals several 
inherent properties of the open system control landscapes. The findings here provide a 
basis to expect that (a) the search for an optimal control in the laboratory may not be 
significantly hindered by the presence of an environment, (b) no traps exist to limit the 
control yield, and (c) in favorable cases, the environment can aid the control outcome. 
These conclusions are consistent with the observed relative ease of attaining control over 
quantum systems in the laboratory, including those of high complexity in condensed 
phase environments. 

A lack of full controllability may turn some of the saddle critical points for control- 
lable systems into local maxima, and thus generate false traps. Moreover, the influence 
of singular extremal controls may also be significant when the system is not controllable, 
because they often appear at the boundary of the reachable set and can become optimal 
solutions to the control problem. Further studies should take into account such limiting 
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Fig 3: Control fields for open quantum control systems: (a) in presence of the environment 
initially at temperature T = 300-fT; (b) in absence of the environment. 

factors as well as (1) uncontrollable degrees of freedom and (2) the employment of re- 
stricted control fields, which might generate apparent false traps in the search for effective 
controls. 
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APPENDIX A: ESTIMATION OF THE NUMBER OF CRITICAL SUBMAN- 
IFOLDS 



An estimate can be given for the number o 



meration is generally NP-hard computable 



33, 



critical submanifolds, whose precise enu- 



38] . The following formula matches well to 



the actual number of critical submanifolds when there are many rows of the contingency 
table and some of the row sums are large. Suppose that r > s, i.e., 6 has higher degen- 



eracies, then we have the Gaussian formula 
of the isolated system: 



39] for the number of critical submanifolds 



S ~ rMexpf-^), (Al) 



where 

Q -2^ E i-—> M = ll DlJ \(s-i)\ ' p = ^ s(s+i) ■ 

j=l i,j 13 v 1 i=l v ' 

Consider the following circumstances: 

(1) For systems with thermal environments at finite nonzero temperature, the marginal 
sums of the contingency tables are 

Dij = d h E k = \e k , i = 1, • • • , r; j = 1, ■ ■ ■ , m; k = 1, ■ ■ ■ , s, 

from which we obtain the approximate number of critical submanifolds associated with 
the A-dimensional thermal environment: 

For large environment dimension A, the number of landscape saddles increases exponen- 
tially with A. 

(2) In the limit of zero temperature, the environmental equilibrium state reduces to a 
pure state. In this case, the number of critical submanifolds remains constant when A is 
sufficiently large, and the precise number is given by the equation ffTTj) . 

20 



(3) In the limit of infinite temperature, = A l I\ becomes a completely mixed state, 
which leads to the marginal sums of the contingency tables: 

A = Adi, Ej = Xej, i = 1, • • • , r; j = 1, • • • , s. 



One can show that for large A, 



_i I s - 1 ) 1 



which leads to 



S2C 2/3 



S'C T— I- 1 



A (r-l)( S -l)_ 



(A3) 



Therefore, the approximate number of critical submanifolds increases polynomially with 
A in the limit of infinite temperature. 
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